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Abstract 


Solutions  of  the  equation  e2u^4^  +  zu(2)  +  otuO)  +  B]U  -  B22zu  =  0,  where 
a,  Bi >  $zZ  are  constants,  e2  <<  1,  and  z  is  the  independent  variable  are  obtained 
using  the  Laplace  integral  technique.  This  equation  describes  the  propagation 
of  high  frequency  electrostatic  waves  near  plasma  resonance  in  a  magnetized 
plasma  with  a  longitudinal  density  gradient  and  is  a  generalization  of  an 
equation  studied  by  Wasow  and  by  Rabenstein  in  the  context  of  boundary  layer 
phenomena.  The  solutions  of  this  fourth  order  equation  in  which  the  associated 
second  order  equation  (i.e.,  e2  =  0)  exhibits  both  a  singularity  (at  z  =  0)  and 
a  turning  point  (at  z  =  B1/S22)  fall  readily  into  two  classes.  One  class  re¬ 
sembles  Airy  functions  and  exists  only  for  e2  not  equal  to  zero.  In  the  other 
class,  the  solutions  are  related  to  confluent  hypergeometric  functions  and  can 
be  viewed  as  solutions  of  the  second  order  equation  with  small  corrections 
proportional  to  e2.  Using  the  integral  representations  of  solutions  it  is 
demonstrated  that  each  class  of  solutions  can  generate  the  other  when  the  inde¬ 
pendent  variable  crosses  the  singular  point.  This  is  the  physical  phenomenon  of 
mode  conversion.  Asymptotic  descriptions  of  both  classes  of  solutions  are  given 
and  the  form  of  the  solutions  near  the  singular  point  is  expressed  as  a  power 
series . 


PACS  number:  52.35.Fp 
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1 •  Introduction 

This  paper  investigates  the  solutions  of  a  fourth  order  differential 
equation  which  arises  in  the  description  of  high  frequency  electrostatic  waves 
near  plasma  resonance  in  a  magnetized  plasma  with  a  zero  order  density  gradient 
along  the  magnetic  field.  The  behavior  of  the  electrostatic  potential  is 
described  by  Poisson's  equation  which  can  be  written  as 

V  •  k  •  V_4>  =  i(co/c)  V_  •  k  •  A  ,  (1.1) 

where  k  is  the  plasma  dielectric  tensor,  $  the  electric  potential,  A  the  vector 
potential,  u>  the  angular  frequency  of  oscillation  and  c  the  speed  of  light. 

A  harmonic  time  dependence  of  the  form  exp(-iwt)  is  assumed.  In  (1.1)  the  term 
containing  the  vector  potential  can  be  viewed  as  a  driving  source  term,  which 
can  be  physically  identified  with  an  externally  launched  electromagnetic  wave, 
as  might  be  the  case  in  a  laboratory  or  ionospheric  experiment.  From  this  point 
of  view  (1.1)  can  be  solved  for  interesting  physical  applications  by  obtaining 
the  appropriate  Green's  function,  a  task  which  requires  knowledge  of  the  solu¬ 
tions  of  the  associated  homogeneous  equation 

V  •  <  •  V<J>  =  o  .  (1.2) 


When  thermal  corrections  associated  with  the  motion  of  plasma  along  the  magnetic 
field  are  retained,  the  plasma  dielectric  tensor  becomes  a  second  order  differential 
operator  and  (1.2)  can  be  written  in  dimensionless  form  as* 
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To  obtain  (1.3)  it  is  assumed  that  the  plasma  has  a  linear  density  gradient 
with  scale  length  L  along  the  magnetic  field  direction.  This  assumption  is 
appropriate  for  many  physical  applications  and  retains  the  important  physical 
processes.  In  (1.3)  u  is  the  electric  potential  suitably  normalized,  and 
z  is  the  distance  along  the  magnetic  field  normalized  to  the  density  scale 
length  L.  The  small  parameter 


e2  =  (kDL)-2  , 


in  which  is  the  Debye  wave  number.  For  typical  ionospheric  plasmas  the 
parameter  e2  can  be  less  than  10-6.  The  other  parameters  in  the  plasma  have 
the  values 


a  =  1  ; 


01  =  (kxL)2fi  2/(a)2  -  a  2)  ; 


e22  =  (w/n  )26!  , 


where  ft  is  the  electron  cyclotron  angular  frequency  and  kx  is  the  fixed  wave 
number  perpendicular  to  the  magnetic  field.  In  obtaining  solutions  of  (1.3)  we 
do  not  restrict  ourselves  to  the  parameter  values  given  in  (1.5) -(1.7).  We  do, 
however,  assume  that  all  parameters  are  real  and  that  t2,  and  g2  are  posi¬ 
tive.  These  assumptions  apply  to  a  plasma  in  which  the  wave  frequency  is 
larger  than  the  electron  gyrofrequency  (i.e.,  w  >  ft_) . 
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Equation  (1.3)  supports  two  distinct  classes  of  solutions:  thermal 
modes  and  cold  modes.  The  first  of  these  classes  represents  short  wavelength 
modes  in  the  sense  that  these  solutions  exist  only  when  e2  (and  hence  k^-2) 
is  not  zero.  The  prototype  equation  for  this  class  is  obtained  from  (1.3) 
by  setting  8^  and  g22  equal  to  zero, 

e2u(1*)  +  zu^2^  +  au^1^  =  0  .  (1.8) 

Equation  (1.8)  approximates  (1.3)  whenever  the  term  e2u^4-*  is  large  in  com¬ 
parison  with  (Biu  -  822zu).  In  this  situation  the  solutions  of  (1.3)  can  be 
obtained  from  the  solutions  of  (1.8)  by  adding  corrections  proportional  to 
Bi  and  2 •  In  the  WKB  sense  the  thermal  class  solutions  are  short  wavelength 
because  the  term  e2u^  in  (1.3)  can  be  large  in  comparison  with  Biu  and 
822zu.  As  shown  in  Sec.  3,  the  solutions  of  (1.8)  are  related  to  Airy  functions 
since  they  are  proportional  to  the  (a  -  2)  derivative  of  Airy  functions  of 
negative  argument  when  a  is  an  integer. 

The  second  class  of  solutions  associated  with  (1.3)  is  comprised  of  cold 
plasma  modes  in  the  sense  that  these  solutions  exist  even  when  z‘  (and  hence 
kp-')  is  equal  to  zero.  The  prototype  equation  for  the  cold  mode  class  is 
obtained  from  (1.3)  by  setting  e2  =  0, 

zu^2^  +  au^^  +  SjU  -  B22zu  =  0  .  (1.9) 

The  solutions  of  (1.9)  are  related  to  confluent  hypergeometric  functions.  When 
the  term  e2u^  is  small,  solutions  of  (1.3)  can  be  obtained  from  solutions  of 
(1.9)  by  adding  small  corrections  proportional  to  e2 .  In  the  WKB  sense  the 


solutions  in  the  cold  mode  class  are  long  wavelength  in  that  c^u'-4-'  is  small 
in  comparison  with  the  other  terms  in  (1.3) 

Solutions  in  the  two  classes  are  generally  distinguished  by  the  disparity 
in  their  wavelengths.  However,  near  plasma  resonance  (z  =  0) ,  the  WKB  wave¬ 
lengths  of  the  two  modes  become  comparable  and  mode  conversion  occurs.  That  is, 
solutions  of  one  class  generate  solutions  of  the  other  class.  This  mode  con¬ 
version  process  is  clearly  exhibited  by  solutions  obtained  in  this  study. 

The  second  order  equation  (1.9)  obtained  from  (1.3)  by  setting  c2  =  0 
exhibits  both  a  singularity  at  z  =  0  and  a  turning  point  at  z  =  3 1 / 62 2 .  The 

existence  of  the  turning  point  distinguishes  (1.3)  from  an  equation  previously 
2  3 

studied  by  Wasow  and  by  Rabenstein.  The  results  obtained  by  these  authors  can 
be  recovered  in  the  limit  $2  *  0*  In  the  plasma  application  the  singularity 
corresponds  to  plasma  resonance  and  the  turning  point  to  upper  hybrid  resonance. 
As  shown  in  Sec.  5  the  existence  of  the  turning  point  profoundly  affects  the 
structure  of  the  cold  plasma  modes  and  allows  for  the  existence  of  solutions 
which  exhibit  no  mode  conversion  for  certain  restricted  parameter  values. 

The  paper  is  organized  as  follows.  In  Sec.  2  we  introduce  integral  repre¬ 
sentations  for  the  solutions  of  (1.3)  and  describe  the  contours  associated  with 
the  solution  set.  In  Sec.  3  the  solutions  corresponding  to  the  thermal  modes 
are  obtained.  In  Sec.  4  expressions  for  the  solutions  corresponding  to  the 
cold  plasma  modes  are  derived.  In  Sec.  5,  the  mode  conversion  process  and  the 
linear  independence  of  the  solutions  is  discussed.  Finally,  in  Sec.  6  the 
principal  findings  are  summarized. 


2 .  General  Properties  of  Solutions 


Since  all  of  the  coefficients  of  u,  and  its  derivatives  in  (1.3)  are 
linear  in  the  independent  variable  z,  general  solutions  may  be  found  in  the 
form  of  Laplace  integrals.  Following  Coddington  and  Levinson,^  but  using  the 
kernel  exp(-sz)  in  the  Laplace  integral  instead  of  exp(sz),  we  obtain  solutions 
of  (1.3)  in  the  form 


where 


s 

u(a,61,02>e2,z)  =  /  (e”sz/P(s) }exp[-/  (Q(t)/P(t) }dt]ds  , 
C 


p(t)  =  t2  -  02 2  ; 


Q(t)  =  e2t4  -  at  +  0X  . 


The  contours  C  in  (2.1)  are  chosen  such  that,  at  the  endpoints  of  the  contours, 
the  following  condition  is  satisfied 

s 

e~SZ  exp[-/  (Q(t)/P(t) )dt]  =  0  . 


Using  (2.2)  and  (2.3)  we  obtain 

s  s  £2t4  -  at  +  0i 

/  (Q/P)dt  =  /  [ - ]  dt 

t2  -  022 

=  e2(s3/3  +  022s)  -  (a/2)2n(s2  -  022) 


3  =  (e2B23  +  8i/62)/2  . 


Inserting  (2.5)  into  (2.1)  the  solutions  can  be  written  as 


where 


u  =  /(s  +  e2)a+  (s  -  B2)a~  exp[-£2s3/3  -  sz]  ds 
C 


z  +  e2B22  , 


a+  =  a/2  -  1  +  g  , 


a  =  a/2  -  1  -  6  . 


The  endpoints  of  the  contours  C  are  chosen  to  satisfy  the  condition 


(s  +  e2)a/2  +  6(s  -  62)a^2  ”  Pexp[-e2s3/3  -  sz]  =  0  • 


It  is  worth  noting  that  the  basic  equation  (1.3)  is  invariant  under  the 
transformation 


(a,  Bj,  62,  e2,  z)  (a,  -Bi,  B2,  -e2,  -z)  . 


Thus  another  family  of  solutions  of  (1.3)  is  given  by 


1  (a,  Bi,  B2,  e2,  z)  =  u(a,  -Bi,  B2,  -e2,  -z) 


=  /,  (s  +  B2)a“(s  -  B2)a+  exp[t2s3/3  +  sz]  ds  , 
C 


where  the  endpoints  of  the  contours  C'  satisfy  the  condition 


(s  +  B2)a//2  6(s  -  B2)a,/2  +  gexp[e2s3/3  +  sz]  =  0 
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In  the  following  we  investigate  the  functions  u  as  defined  by  (2.7)  in  detail 
and  use  the  transformation  (2.12)  to  obtain  the  functions  <f>. 

The  condition  (2.11)  that  must  be  met  at  the  endpoints  of  the  contours 
can  be  satisfied  at  large  s  by  choosing  the  real  part  of  e2s3/3  to  be  positive. 

Since  we  have  chosen  e2  real,  this  requirement  becomes,  with  j  =  1,2,3, 

2tt (2- j ) / 3  -  tt / 6  <  arg  s  <  tt/6  +  2tt (2- j )/3  .  (2.15) 


Thus  as  shown  in  Figure  1,  there  are  three  open  sectors  i  s-plane  of 

angular  width  tt/ 3  centered  about  2tt/3,  0,  and  -2ir/3  in  which  the  contour  may 
go  to  infinity  and  satisfy  the  condition  (2.11).  Since  we  have  stipulated  that 
a,  Bj,  82  and  £2  are  real  and  positive,  the  quantity  a/2  +  8  is  real  and  positive. 
Thus  the  condition  (2.11)  can  also  be  satisfied  at  the  point  s  =  -82  because 
the  quantity  (s  +  82)a/,2+^  is  then  always  zero.  Since  the  quantities  a+  and  a_, 
while  real,  are  not  necessarily  integers,  the  points  s  =  ±02  are,  in  general, 
branch  points.  If  the  s-plane  is  cut  from  82  to  positive  infinity  and  from 

-82  to  negative  infinity  along  the  real  axis  as  shown  in  Figure  1,  the  multi¬ 
valued  function  (s  -  82)a~(s  +  32)M+  in  the  integrand  of  (2.7)  can  be  written  as 

,  „  ,a_,  „  ',«+  r  i6i>a_,  i0?,a+ 

(s  -  62)  (s  +  62)  +  =  (ri  e  x)  (r2  e  2) 


=  exp[a_2,nri  +  ia_(0j  +  2nn)  +  a+£nr2  +  ia+(02  +  2irm)]  , 


(2.16) 


with  - 7t  <  0 2  <  tt  and  0  <  ©i  <  2tt .  In  (2.16)  the  various  branches  of  the  £n 
functions  are  represented  by  the  integers  m  and  n  which  denote  various  sheets 
of  the  Riemann  surface.  Each  sheet  of  the  Riemann  surface  can  thus  be  labeled  by 
the  pair  of  integers  (m,n).  The  s-plane  represented  in  Figure  1  corresponds  to 
the  principal  branches  of  the  £n  functions  m=n=0  and  is  labeled  by  (0,0). 
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The  endpoints  for  contours  corresponding  to  solutions  of  (1.3)  must  he 
chosen  to  satisfy  the  relation  given  in  (2.11).  This  condition  can  be  met  by 
choosing  the  contours  that  begin  and  end  at  infinity  in  any  of  the  sectors 
j  =  1,2,3  or  contours  that  begin  at  the  branch  point  s  =  -f32  ar>d  end  at 
s  =  -32  or,  alternatively,  proceed  to  infinity  within  the  numbered  sectors. 
Contours  which  correspond  to  solutions  in  the  thermal  mode  class  are  shown  in 
Figure  2.  These  contours  begin  at  infinity  in  one  sector  and  end  at  infinity 
in  another  sector.  The  solutions  obtained  from  (2.7)  by  integrating  along 
these  contours  are  labeled  where  j  refers  to  the  sector  opposite  the  contour. 
For  example,  the  contour  beginning  in  sector  2  and  ending  in  sector  1  corres¬ 
ponds  to  the  solution  labeled  A3.  The  contour  for  the  solution  A2  begins  in 
sector  1  and  crosses  the  branch  cut  along  the  negative  real  axis  before  pro¬ 
ceeding  to  infinity  in  sector  3.  Tlws  part  of  the  contour  for  the  solution  A2 
lies  in  the  (1,0)  Riemann  sheet  (i.e.,  the  sheet  with  m  =  1,  n  =  0)  and  is  thus 
shown  as  a  broken  line  in  Fig.  2. 

Contours  corresponding  to  solutions  in  the  cold  plasma  mode  class  have  at 
least  one  end  point  at  s  =  -82  and  are  shown  in  Figs.  3(a)  and  3(b).  The  contour 
corresponding  to  the  solution  labeled  B0  starts  at  the  point  s  =  -B2,  encircles 
the  point  s  =  82  in  the  counter-clockwise  direction  and  ends  again  at  s  =  -B2* 

The  contour  for  B0  crosses  the  branch  line  along  the  positive  real  axis  and 
thus  crosses  onto  the  (0,1)  Riemann  sheet  (i.e.,  m  =  0,  n=  1).  The  portion  of 
the  contour  lying  in  the  (0,1)  sheet  is  shown  as  a  dashed  line.  Contours 


corresponding  to  solutions  labeled  B^  where  j  corresponds  to  the  sector  in 
which  they  proceed  to  infinity  are  illustrated  in  Figure  3(a).  These  contours 


start  at  the  point  s  =  -f?2  and  proceed  to  infinity  in  sector  j  in  such  a  fashion 


that  the  radius  vector  from  s  =  82  to  a  point  on  the  contour  moves  in  the 
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clockwisc  direction  as  the  point  proceeds  to  infinity.  Thus  the  contour  for  the 
solution  B;>  passes  above  the  point  s  =  3 2-  Contours  for  solutions  labeled  11. 
are  shown  in  Figure  3(b).  These  contours  also  start  at  s  =  -3;>  but  the  radius 
vector  from  s  =  -t>2  to  a  point  on  the  contour  moves  in  the  counter-clockwise 
direction  as  the  point  proceeds  to  infinity  in  sector  j.  Thus  the  contour  for 
B2  passes  below  the  point  s  =  32.  Note  that  the  contour  for  B]  shown  in 
Figure  3(b)  crosses  the  branch  cut  along  the  positive  real  axis  and  thus  passes 
onto  the  (0,1)  Riemann  sheet.  The  portion  of  the  contour  on  the  (0,1)  sheet 
is  again  shown  as  a  dashed  line.  While  not  illustrated  in  Figure  3(a),  the 
contour  for  the  solution  B3  would  similarly  cross  the  branch  cut  along  the 
positive  real  axis  but  in  a  clockwise  direction  and  thus  would  pass  onto  the 
(0,-1)  Riemann  sheet.  The  contours  for  the  solutions  B^  and  B^  (j  =  1,2,3) 
proceed  to  infinity  along  the  same  asymptotic  direction  in  sector  j . 

Although  illustrated  only  for  the  principal  sheet,  the  contours  corresponding 
to  solutions  Bj ,  B ^ ,  BQ  and  may  begin  on  any  sheet  of  the  Riemann  surface. 

In  order  to  distinguish  among  functions  corresponding  to  contours  on  different 
sheets  of  the  Riemann  surface  we  introduce  the  notation  X(m,n;p),  where  X  is 
any  of  the  solutions  A^ ,  B0 ,  ,  B^ ,  m  and  n  are  integers  specifying  the  sheet, 

and  p  denotes  dependence  on  the  parameters  a,  3j ,  32 »  e2,  z.  In  this  notation 
the  values  m  and  n  indicate  the  sheet  on  which  the  contour  begins.  As  described 
above  some  contours  begin  on  one  sheet  and  end  on  an  adjacent  sheet.  Thus 
A?(0,();p)  indicates  the  contour  beginning  in  the  sheet  with  m=0,  n=0  but,  as 
shown  in  Fig.  2,  it  ends  in  the  sheet  with  m=l,  n=0.  For  brevity  of  notation  we 
omit  the  dependence  on  m  and  n  or  other  parameters  unless  they  are  needed  to 
clarify  the  discussion.  Furthermore,  if  the  values  of  m  and  n  are  not  explicitly 
indicated  the  principal  values  m=0,  n=0  are  to  be  assumed.  As  an  illustration 


we  have  drawn  in  Figure  3(a)  the  contour  corresponding  to  the  function  B3(0,1). 
The  contour  begins  on  the  sheet  (0,1)  (and  is  thus  shown  dashed)  and  passes 
onto  the  principal  sheet  when  it  crosses  the  branch  cut  extending  from  82 
positive  infinity.  Finally,  functions  corresponding  to  contours  on  different 
sheets  are  simply  related.  Employing  the  notation  described  above  and  using 
(2.16)  one  can  write 


v ,  .  i27rma.  i27rna_ 

X(m,n)  =  e  +  e  X(0,0) 


(2.17) 


Since  the  integrand  in  (2.7)  is  analytic  throughout  the  entire  s-plane 
except  at  the  points  s  =  ±8?,  Cauchy’s  theorem  can  be  used  to  establish  relation¬ 
ships  among  the  various  solutions  A^  ,  B0,  B^ ,  and  B^  by  using  combinations  of 
the  appropriate  contours.  In  order  to  facilitate  the  derivation  of  relations 
among  the  solutions  all  of  the  contours  in  Figs.  2  and  3  have  been  combined  into 
Fig.  4.  Refering  to  Fig.  4  and  using  Cauchy's  theorem  the  following  relationships 
can  be  established 


B2  =  Bj  -  A 3 


63  =  B2  -  A3 (0 , 1) 


Bi  =  B3(1,0)  -  a2 


B3(0,1)  =  B2(0,1)  -  Ax 


B3  =  B2  -  Aj 


(2.18) 


(2.19) 


(2.20) 


(2.21) 


(2.22) 


B.  -  B. (0,1)  =  B  ;  (j  =  1,2,3) 


(2.23) 
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The  relations  (2.18),  (2.20)  and  (2.22)  can  be  readily  verified  by  referring 
to  Fig.  4.  To  verify  the  remaining  relations  it  is  helpful  to  picture  the 
contours  for  §2  or  B2(0,l)  as  crossing  the  branch  cut  which  extends  from  82 
to  infinity.  The  relations  (2. 18)- (2. 23)  prove  useful  in  obtaining,  among 
other  things,  analytic  continuations  of  the  functions  and  B^ .  Having 
determined  the  set  of  contours  which  yield  solutions  of  (1.3)  in  the  integral 
representation  (2.7)  we  next  proceed  to  evaluate  these  solutions  in  detail. 

We  first  examine  the  solutions  A.  before  proceeding  to  investigate  the  solutions 


B„ ,  B .  and  B . . 
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3.  The  Solutions  A. 

_ I 

The  functions  A^(z)  are  defined  by  the  integral  representation  (2.7)  in 
which  the  paths  of  integration  are  given  by  the  contours  C(A^.)  as  illustrated 
in  Fig.  2.  In  this  section  we  first  discuss  the  asymptotic  behavior 
of  the  functions  A_.  (z)  for  |z|  »  e‘'^3  and  then  derive  a  power  series  expansion  which 
is  particularly  useful  for  small  argument,  z  -»■  0.  In  both  cases  it  is 
convenient  to  introduce  the  transformation 


a  =  e2/3 s;  n  =  e~2/3z;  oq  =  0(^2)  =  e2//>3S2 


(3.1) 


where  z  =  z  +  e2022  ~  z  when  e2  <<  1.  Making  use  of  (3.1)  in  (2.7)  we  obtain 
the  integral  representation 


Aj  "»  ■  E 


■  2  (a-1  )/3 


/  (a  -  a  )a-  (a  +  a  )a+  exp[-(a3/3  +  an)]  da  ,  (3.2) 


CCA.) 


where  C'(Aj)  is  the  image  of  the  contour  C(A^)  under  the  transformation  (3.1) 

To  evaluate  (3.2)  asymptotically  we  employ  the  saddle  point  method  of 
integration.  We  note  that  the  derivative  of  the  exponent  f(a)  =  -(a3/3  +  an), 
namely  /'(a)  =  -(a2  +  n),  when  equated  to  zero  yields  two  roots,  or  saddle  points, 

t, 

at  a^  =  +in 2.  For  our  purpose  we  discard  the  minus  sign  and  choose  the  saddle 
point 


.  \  1  1  i0  11  . 

0^=10;  n  =  | n I e  ;  |n|  >>  1  , 


(3.3) 


■ 
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from  which  we  deduce  that,  when  0  =  0,  =  in2  with  n  real  and  positive;  when 

9  =  TT>  ag  =  eilT  Ini’5  =  -(-n)”5  with  n  real  and  negative,  and  finally,  when 
9  =  2tt,  os  =  e37T^2  |n|^  =  -in^,  with  n  again  real  and  positive.  These  values 
of  represent  the  principal  saddle  points  and  their  location  in  the  cr-plane 
is  illustrated  in  Fig.  5.  The  principal  saddle  points  at  0  =  0,  tt  and  2ir 
correspond  respectively  to  the  contours  for  A3,  A2 ,  and  A3.  The  other  choice 
of  sign,  =  -in  ,  yields  nothing  new. 

The  contours  C'(A^)  in  (3.2),  as  shown  in  Fig.  5,  are  asymptotic  to  the  rays 
with  phase  0,  2tt/3,  and  4ir/3.  Accordingly,  making  use  of  (3.3)  one  has  the 
following  argument  ranges 

2tt  (2- j  )/3  +  2ir/3  <  arg  a g  <  4ir/3  +  2tt  (2- j) /3 ; 

4ir(2-j)/3  +  tt/3  <  0  <  5ir/3  +  4ir(2-j)/3  .  (3.4) 

As  shown  in  Fig.  5,  we  can  use  (3.4)  to  trace  the  path  of  the  saddle  point  with 
|  n  |  fixed  as  0  varies  from  -it  to  3tt,  or  arg  a  varies  from  0  to  2ir.  It  is  seen 
from  Fig.  5  that  the  chosen  saddle  point  in  (3.3)  corresponds  to  the  func¬ 

tions  A3,  A2,  and  Aj  and  traces  a  circle  of  radius  0  =  Ini’5.  However, 
we  note  that  starts  just  above  the  right  hand  branch  cut  of  Fig.  1  and  traces 
a  semi-circle  in  the  upper  half  of  the  a-plane,  at  which  point  a crosses  the 
left-hand  branch  cut  from  the  principal  sheet  onto  the  adjacent  sheet  of  the 
Riemann  surface,  m  =  1 ,  n  =  0  in  accordance  with  (2.16),  and  traces  a  semi-circle 
on  the  lower  half  of  the  (1,0) sheet. 

We  observe  from  (3.4)  that  the  path  of  steepest  descents  for  A3  lies  entirely 
on  the  principal  sheet  of  the  Riemann  surface  and  therefore  the  saddle  point 
integration  of  (3.2)  is  independent  of  the  branch  cuts,  which  is  to  say 
that  the  multivalued  factors  (o  -  °0)a"  and  (o  +  °Q)a+  in  the  integrand  of  (3.2) 
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are  assigned  their  principal  values.  For  A2  the  path  of  steepest  descents 
starts  at  infinity  on  the  principal  sheet  and  terminates  at  infinity  on  the 
adjacent  m  *  1,  n  =  0  sheet  thus  crossing  the  left-hand  branch  cut.  However, 
since  the  saddle  point  also  crosses  onto  the  (1,0)  sheet,  the  saddle  point 

integration  is  carried  out  assuming  that  the  multivalued  factors  in  (3.2) 

vary  continuously  as  the  saddle  point  crosses  the  branch  cut.  Finally, 
as  regards  the  path  of  steepest  descents  for  Aj ,  we  note  that  it  is  also 
independent  of  the  branch  cuts  but  now  lies  entirely  on  the  m  =  1,  n  =  0 
sheet  and,  therefore,  the  saddle  point  integration  of  (3.2)  evaluates  Ai(l,0). 

The  multivalued  factor  is  evaluated  according  to  (2.16)  with  m  =  1,  n  =  0. 

Next,  we  wish  to  determine  the  direction  of  traversal  through  the  chosen 
saddle  point  a  .  In  the  vicinity  of  a  we  can  write  -x2  =  %  /"(a  ) (a  -  o  )2  +  .... 
where  x  is  real  and  positive  (x  >  0)  after  passing  through  the  saddle  point. 

In  the  present  instance  -H  f"  (og)  ~  °s  ~  ih2.  Hence,  extracting  the  square 
root  and  putting  w  =  a  -  a  we  have  x=  ±l-h  f"  (o  )]**  w  =  +(o  )**  w.  Ignoring 

the  plus  sign,  writing  -1  =  e"17T  and  putting  arg  x  =  0  yields 

arg  w  =  -rr  -  arg  (a  )  2  =  3n/4  -  0/4  ,  c 

which  says  that,  for  A3  and  0  =  0,  arg  w  =  3tt/4 ;  for  A2  and  0  =  it,  arg  w  =  tt/2, 
and  for  Ax  and  0  =  2tt,  arg  w  =  ir/4.  We  note  that  (3.5)  gives,  for  A3  and  At , 
the  correct  direction  of  traversal  through  the  saddle  point  in  accordance  with 
the  arrows  drawn  in  Figs.  2  and  5.  However,  for  A2  the  direction  of  traversal  is 
opposite  to  the  direction  of  the  arrow,  which  gives  a  minus  sign  to  be  attached 
to  the  final  result  for  A2 . 


Asymptotically,  |z|  »  e  7  and  with  z1  <<  1  yields  z  =  z  +  e2$2  -  z 
-2/3 

and  n  -  e  z.  Furthermore,  from  (2.6)  8  =  e2823/2  +  81/282  =  81/282  =  80, 
from  which  instead  of  the  exponents  a+  and  a_  defined  by  (2.9)  and  (2.10) 
we  introduce 


a+  =  a/2  -  1  +  80  ;  a_  =  a/2  -  1  -  80  . 


After  this  preamble,  we  can  now  apply  the  familiar  leading  term  formula  of  the 
method  of  steepest  descents  to  compute  an  asymptotic  representation  for  the 

functions  (n)  given  by  the  integral  representations  (3.2).  Thus  we  obtain, 

v 

for  the  saddle  point  =  in2  and  j  =  1,2,3,  the  leading  terms: 


Aj(y,0;n)  ^  /iT  (-)je_2(:a"l)/3(as  -  0Q)a-(as  +  gQ)a+exp[f (ag)]  ^  (1+0(0,  3)> 


Jts  (-)^e  2^a  1^3(in^-cr  )a-(in^+a  )a+  e~a^  (1+0(1/C)> 


[in^r 


where 


C  s  (2/3)n~ 


(2/3)e‘V/2  ;  Y  =  -  (-)J]Y  , 


0  Im(os)  >  0 

1  Im(o  )  <  0 


We  have  introduced  the  factor  y  to  insure  the  continuity  of  the  function  A2  as 
the  saddle  point  crosses  the  branch  cut  along  the  negative  real  axis.  The  general 


solution  Aj (m,n;n)  can  be  found  from  (3.7)  by  using  (2.17). 

From  (3.8)  we  deduce  that  for  j  =  3  and  0  =  0,  -i£  =  -(2/3)in3//  with  n 

3/2 

real  and  positive;  for  j  =  2  and  6  =  tt ,  -ic  =  -(2/3)(-n)  '  with  n  real  and 

.  3/2 

negative,  and  finally,  for  j  =  1  and  8  =  2tt,  -i c  =  (2/3) in  with  n  again  real 

and  positive.  Finally,  to  conclude  our  asymptotic  leading  term  presentation  we 

2/3 

translate  (3.7)  into  a  function  of  z.  Recalling  that  oq  =  e  g2,  we  obtain, 
after  some  algebraic  manipulations. 


Aj  (Y,0;z)  =  (~)J  (iz 2/e-B?)a~ (iz^/e+g2)a'1'  Y  U+0(l/rJ}  , 


+  /rT 


■U 


[iez^] 


(3 


where  the  multivalued  factor  is  evaluated  as  discussed  above. 

In  (3.7)  we  have  retained  the  term  eleven  though  we  have  assumed  |os  |  >>  oq 
(i.e.,  |z|  >>  e2822)  in  order  to  preserve  the  topological  structure  of  the  cut 
o-plane.  It  is  useful  to  obtain  a  form  for  A^  in  which  oq  is  ignored  in  comparison 
with  but  this  procedure  will  alter  the  topology  of  the  cut  o-plane.  Ignoring 
the  oq  term  in  (3.7)  is  formally  equivalent  to  setting  oq  =  0  so  that  there  is  a 
single  branch  point  at  the  origin  rather  than  two  branch  points.  The  multi- 

ct  ci  ct-  2 

valued  factor  in  (3.7)  becomes  og  "o  +  =  °s  •  Taking  the  branch  cut  to  lie 
along  the  positive  real  axis  and  the  principal  sheet  to  have  argument  range 
0  <  arg  og  <  2rr ,  we  see  that  the  argument  range  for  the  multivalued  factor 
(o  -  o  )  “(a  +  a  )a+  in  (3.7)  exactly  coincides  with  that  of  a  a"a  a+  when  a  =0. 

Noting  from  (2.17)  that  Ai(l,0;n)  =  e1?Tra+A](n)  we  can  then  write 


Aj(n)  ^  (-)■ 


A 


r2(a-l)/3(in2)a-5/Vif' 


where  the  phase  factors  are: 


P3  =  P2  =  1  ;  Pi 


-i2ira+ 

e  + 


(3.11) 


( 


In  (3.11)  the  argument  range  of  (in  )  is  from  0  to  2tt . 

As  mentioned  previously  the  functions  (n)  are  closely  related  to  Airy 
functions.  This  relationship  is  best  elucidated  if  we  consider  the  case  =  0, 

which  is  equivalent  to  setting  82  =  81  =  0,  with  the  result  that  (1.3)  reduces 

2/3 

to  (1.8).  Putting  n  =  e  z  in  accordance  with  (3.1)  and  writing  x(n)  instead 
of  u(z),  we  obtain  from  (1.8)  the  equation  x  nx  ^  +  (a  -  l)x  =  0.  When  a 
is  an  integer,  it  can  be  shown  that  x(n)  is  proportional  to  the  (a  -  2)  derivative 
of  any  Airy  function  of  negative  argument,  i.e.,  any  solution  of  vv  '+  ryv  =  0. 

When  (o  -  2)  is  a  negative  integer  the  above  statement  must  be  interpreted  as 
the  | a  -  2 |  integral  of  v(n).  The  same  conclusions  follow  from  the  integral 
representation  (3.2)  upon  setting  oq  =  0;  that  is,  A^ (n)  in  (3.11),  when  a  is 
an  integer,  is  proportional  to  the  (a  -  2)  derivative  of  the  asymptotic  leading 
term  of  an  Airy  function  of  negative  argument,  and  hence  is  a  solution  of  our 
basic  equation  (1.3). 

We  now  proceed  to  discuss  the  small  argument  behavior  of  the  function  A^(n). 
For  this  purpose  we  make  use  again  of  the  scale  transformation  (3.1)  and  the 
integral  representation  (3.2)  wherein  a+  and  a  have  their  original  definitions 
(2.9)  and  (2.10).  To  evaluate  (3.2)  we  first  construct  the  power  series  expansion 


(a  -  a  )  (0+0) 


+  _ 


ct-2  r- 

°  l 


m 


c  w 


1 

1 


wl  <  1  . 


( 
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where  w  =  aQ/v  and  the  expansion  coefficients 


c 

m 


co 


l 

n=0 


- 

ot. 

<*+ 

n 

m-n 

J 

l 

(3.14) 


are  given  in  terms  of  the  familiar  binomial  coefficients.  The  expansion  on  the 
right  of  (3.13)  is  absolutely  convergent  for  |wj<l  because  it  arises  from  the  multi¬ 
plication  of  the  two  absolutely  convergent  binomial  expansions  corresponding  to  the 
factors  (a-o^)  ~  and  (o +aQ)  +  .  This  result  follows  from  Cauchy's  Theorem  on  the 
multiplication  of  absolutely  convergent  series. ^  Since  the  power  series  (3.13) 
is  absolutely  convergent  for  |w|  =  l-<5,  where  6  is  an  arbitrarily  small  number,  it 
also  follows  that  the  series  in  question  is  uniformly  convergent7  for  |w|  <  1  -  <5. 
Hence,  we  can  replace  the  product  (a  -  o^)  ~  (a  +  ctq)  +  in  the  integrand  of  (3.2) 
with  the  uniformly  convergent  series  (3.13)  and  integrate  term  by  term  to  obtain, 
for  j  =  1,2,3, 


Aj(n) 


(a- 1 ) / 3  r 
m=0 


gj  (n,a-m) 


(3.15) 


where  the  phase  factors  Pj  are  given  by  (3.12),  and  where 


gi(n,a)  = 

J  j 

C'(A.) 


,  a-2  -a3/3-an 

da  a  e 


(3.16) 


are  the  same  functions  introduced  by  Rabenstein.  In  (3.16)  the  contour  C'(A^) 
is  chosen  such  that  |oj  >  oq  everywhere  along  the  contour,  which  is  the  condition 


for  the  convergence  of  the  series  (3.13). 
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A  convenient  expression  for  the  functions  gj(n,ct),  j  =  1,2,3,  is  obtained 
by  expanding  the  factor  exp[-an]  into  a  power  series 


-an  v  r  r  -.k „  . n  k  k 
e  =  I  [(-)  /k!]a  n 

k=0 


(3.17) 


which  is  uniformly  convergent  in  the  finite  a  plane,  |a|  <  Introducing 
(3.17)  into  (3.16)  and  integrating  term  by  term  we  obtain 


gj (n,a) 


-  I 


[(-)k/k!]nk 


k=0 


,  a-2+k  -a3/ 3 

da  a  e 


C'(A.) 


(3.18) 


Putting  y  =  a3/ 3,  the  integrals  in  (3.18)  become 


I.(a,k) 


da  a“-2+k  e‘a3/3  =  3(a+k"4)/3 


dy  e"X  y(“+k-4i/3 


C'(A..) 


H. 

3 


(3.19) 


where  the  contour  1L  is  the  image  of  C ’  (A_. }  under  the  transformation. 

To  ascertain  the  shape  of  H.  consider,  for  example,  C'(A3)  which  begins  at 
infinity  with  arg  a  =  0  and  proceeds  to  infinity  with  arg  a  =  2tt/3.  Under  the 
transformation,  y  =  a3/ 3,  the  contour  H3  begins  at  infinity  with  arg  y  =  0, 
encircles  the  origin  in  the  counter-clockwise  direction,  and  proceeds  to  infinity 
with  arg  y  =  2n.  Therefore,  H3  is  simply  Hankel's  contour  and  we  can  write 


1 3 (a,k)  =  3(a+k'4)/3 


(0-0 


dy  e~y  y(a+k-4)/3 


(3.20) 


which  can  be  evaluated  in  terms  of  gamma  functions. 
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Proceeding  in  an  entirely  analogous  fashion,  we  ascertain  that  the  contour 
H2  starts  at  infinity  with  arg  y  =  2v ,  encircles  the  origin,  and  terminates  at 
infinity  with  arg  y  =  4ir ,  whereas  the  contour  Hj  begins  with  arg  v  =  4tt  and  ends 
with  arg  y  =  6u .  Thus,  introducing  the  abbreviation 

P  =  p(a,k)  =  (4-ct-k)/3  ,  C3.21) 

9 

and  making  use  of  a  generalization  of  Hankel's  integral  representation, 

Ij  (ot , k)  =  27Ti  expUiCl  -  p)(7  -  2j)]/{3pr(p) }  .  (5.22] 

Making  use  of  (3.18)  and  (3.19)  we  obtain  the  expansion 

CD 

g,(n,a)  =  l  [(-)k/k!]l  (a,k)nR  ,  (3.23) 

3  k=0  3 

where  the  functions  I^(a,k)  are  given  by  (3.22).  Finally,  using  (3.15)  we 
obtain  the  power  series  expansions 

A  (n,a)  =  p  e'2(a'1)/3  l  l  c  a  m[(-)k/k!]l  (a  -  m,k)nk  ,  (3.24) 

3  3  m=0  k=0  3 

where  the  p ^  are  the  phase  factors  given  in  (3.12)  These  functions,  for  j  =  1,2,3, 
correspond  to  their  asymptotic  counterparts  (n)  given  by  (3.7)  and  are  subject 
to  the  same  argument  ranges  for  n  as  given  by  (3.4). 
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4.  The  Solutions  B.,  B.,  and  Bn 
_ J  J  ° 

The  functions  B^  ,  B^ ,  and  B0  are  solutions  of  (1.3)  and  are  given  by  the 
integral  representation  (2.7),  where  the  corresponding  paths  of  integration  are 
shown  in  Figs.  3(a)  or  3(b).  Since  we  are  unable  to  directly  evaluate  these 
integrals,  we  resort  to  a  perturbation  expansion  based  on  the  fact  that  e2  <<  1. 

Thus  expanding  (2.7)  in  a  power  series  in  e2  we  show  in  the  Appendix  that,  to 
first  order  in  e2  , 

ufa.Bi ,B2,e2;z)  =  u(a , B1 , 62 ,0 ; z)  +  e2(3  2u) 

e2=0 

=  uQ  -  e2[h(uQ)/462]  ,  (4.1) 

where  u(a,Bi  ,62,e2 ’,t)  is  any  of  the  solutions  BQ,  or  B^  and  uQ  is  the 
corresponding  solution  evaluated  at  e2  =  0.  The  expression  for  h(u0)  is  given 
in  (A. 2. 6)  in  terms  of  uQ  evaluated  at  shifted  values  of  the  parameters  a  and  Bi . 

We  observe  that  the  functions  uQ (z)  are  solutions  of  the  reduced  (e2  =  0) 
second  order  differential  equation  (1.9)  and  are  given  by  the  integral  represen¬ 
tation 
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xg*-2-*  +  (a  -  x)g^  -  (a/2  -  B0)g(x)  =  0  ,  (4.3) 

which  has  two  independent  solutions  that,  for  our  purpose,  we  write  as 

g^(x)  =  M(a/2  -  B0,a;x)  ;  (4.4) 

g?(x)  =  U (a/2  -  B0,a;x)  ,  (4.5) 

where  M  and  U  are  Kummer  functions.  Since  the  properties  of  these  functions 
are  well-known,  we  proceed  to  obtain  expressions  for  ,  b_.  ,  and  b0  in  terms 
of  Kummer  functions. 

The  function  b0  obtained  from  (4.2)  by  integrating  along  the  contour 
associated  with  B0  is  related  to  the  Kummer  function  M(b,c;z).  This  can  be 
demonstrated  by  using  the  transformation 

s  =  202 (t  -  h)  (4.6) 

in  (4.2)  to  obtain 

b0  =  (232)“' 1  e^2  /  ^  (t  -  I)3"  ta+  e"2B2Zt  dt  .  (4.7) 

0 

The  contour  in  (4.7)  is  the  image  of  C(Bq)  under  the  transformation  (4.6).  It 
starts  at  the  origin,  circles  the  point  t  =  1  in  the  positive  direction  and 
returns  to  the  origin.  The  integrand  in  (4.7)  is  bounded  everywhere  along  the 
contour  for  jz|  <  00  so  that  an  expression  can  be  obtained  for  b  that  is  valid 
for  the  entire  z-plane,  as  is  shown  in  the  following. 
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The  integral  representation  for  the  Rummer  function  M(b,c;z)  for 


Re(b)  >  0  is' 


M(b  c  z)  =  TJsUS^zSin. 

27Ti  r(b) 


Zt  4-b_1  /> 

e  t  (t 


l)C_b_1  dt 


(4.8) 


From  (3.6)  one  has  a+  =  a/2  +  80  -  1,  where  80  =  81/282  is  8  evaluated  at  e2  =  0, 
and  since  a,  81,  82  are  taken  as  real  and  positive,  Re(a/2  +  80)  >  0  so  that 
(4.8)  can  be  used  in  (4.7)  to  obtain 


b°  ■  eC/2  |  J 


(4.9) 


where  a  =  a/2  -  80  and  £  =  282Z.  Finally  using  the  Rummer  transformation 


yields 


M(b,c;z)  =  e  M(c-b,c‘,-z) 


U  _  I01*1  _~C/ j  2-rri  T  (a  -  a)  M,_  .  ^  > 

bo  -  2^2  e  i  rl^T)— (—--ITM(a’a;C) J 


(4.10) 


(4.11) 


An  expression  for  BQ ,  accurate  to  order  e2,  can  thus  be  obtained  by  combining 
(4. 1),  (4. 11)  and  (A. 2. 6)  : 


=  C2S2 )“•'  e<'2  [[1  - 


(e2/ 1202  )  [M(a  +  l,a  +  4;f.)  -  M(a  +  3, a  +  4 ;  C )  ] 


(e2/ 1 2 32  +  c282/4)[M(a  +  2, a  +  2;£)  -  M(a,a  +  2;fj]  1 


(4.12) 


From  (4.12),  the  value  of  B0  can  be  found  to  order  c 2  for  any  z  with  |z|  <  «=. 
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The  solutions  B.  and  B.  are  related  to  a  linear  combination  of  the  Kummer 
J  J 

functions  N1(b,c;z)  and  U(b,c;z).  The  contours  associated  with  the  solutions 
B_.  and  B^  go  to  infinity  in  the  sector  j  in  accordance  with  (2.15)  and  the 

Ct  —  2  —  Z  S 

integrand  of  u0  as  given  in  (4.2)  has  the  asymptotic  value  s  '  e  .In  order 
to  insure  that  the  integrand  remains  bounded  for  all  a  one  must  require  that 
| z |  >  | z0 |  >  0  (for  some  given  z  ) ,  and 


- tt/ 2  <  arg(sz)  <  tt/2 


(4.13) 


However,  the  subscript  j  already  implies  an  argument  range  for  the  variable  of 
integration  s  along  the  contours  associated  with  the  solutions  B^  and  B^ .  As 
given  by  (2.15)  these  ranges  are,  for  the  solutions  B^ , 


2tt  (2  -  j  )/3  -  tt/6  <  arg  s  <  tt/6  +  2ir(2  -  j)/3  . 


(4.14) 


From  (4.13)  and  (4.14)  one  finds  that  the  solutions  b^  are  well-defined  only 
over  certain  regions  of  the  z-plane,  namely. 


2ii  (j  -  2)/3  -  2it/3  <  arg  z  <  2tt/3  +  2n  (j  -  2 ) / 3  . 


(4.15) 


Using  (4.2)  and  the  transformation 


s  =  2B2(t  ♦  h)  , 


(4.10) 


yields 


iO 

°°G 

b.  =  (2e2)a_1  e</2  f  ta_(t  +  l)a+  e'St  dt  , 
J  -1 


(4.17) 


where  the  angle  0  lies  in  the  sector  j.  The  paths  of  integration  given  in  (4.17) 


are  the  image  of  the  contours  corresponding  to  the  solutions  B^  under  the  trans¬ 
formation  (4.16).  Similarly,  from  the  definition  of  the  contours,  we  have 


b.  =  (2  B2 )  °t_  1  e</2 


I 

-1 


»e 


i(0  +  2n) 


ta-(t  +  l)a+  e‘Ct  dt  . 


The  integral  representation  for  the  Kummer  function  U  is 


.  12 


U(b,c;y)  =  (1/2-n-i) r (1  -  b)e 


-bfri 


(0+) 

ie 


tb-‘ci  *  t)c-b-'  e-yt  dt 


°°e 


In  the  integral  of  (4.19)  the  contour  starts  at  infinity  with  argument  0,  circles 
the  origin  in  the  positive  direction  (i.e.,  cuts  the  negative  t-axis)  and  returns 
to  infinity  with  argument  0  +  2ir.  Since  the  integrand  of  (4.19)  is  analytic  in 
the  region  between  the  origin  and  t  =  -1  the  contour  can  be  extended  to  include 
the  point  t  =  -1;  hence 


(0+) 

'  i0 
o°e 


-1 


; 

me 


i0 


i  (0  +  2tt) 
«>e 

. 

-1 


Using  (4.20)  and  (4. 17) - (4 . 19)  yields 


b. 

J 


(262) 


a-l  2irie 


lira 


r Cl  -  a) 


-E/2 

e  U  U(a,cx;0 


Using  (2.23)  together  with  (2.17)  one  finds  that 


B.  =  [ei27Ta  -  l]'1 
.1 


(B.  -  B.)  -  B01 
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which,  for  e-  =  0,  becomes 


b.  =  [ei2lTa  -  l]"1  (b.  -  b.)  -  b0  \ 


Using  (4.21)  to  replace  the  combination  b^  -  b^  in  (4.23)  and  employing  the 


identity  F(a)r(l-a)  =  n/sin  (an)  gives 


.  -nra 

b.  =  (262)a_1r(a)e  ’/2  U(a,a;C)  -  - 2 — 

-1  2i  sin(na) 


The  relation  (2.23)  for  e2  =  0  gives 


I2na 


bj  =  b^  (0,1)  +  b0  =  e  bj  +  b0  . 


Using  (4.25)  to  replace  b^  by  b^  in  (4.24)  gives 


„  /  -nra  , 

,  ^a-1  i2ira„,  ..  -C/2  ^  e  b0 

b.  =  (2B2)  e  r (a)e  U(a,a;C)  -  y 


2i  sin(na) 


With  (4.11)  the  relations  (4.24)  and  (4.26)  can  be  used  to  express  K  and  y 


entirely  in  terms  of  the  Kummer  functions  U  and  M, 


bj  (C) 


=  (202) 


nnd 


(4 .23) 


(4.24) 


(4.25) 


(4.26) 


1 


Vf) 


(2B2)C‘"1  e'^2{r (a)exp (2iria)U(a ,a ; C)  -  si^f^a)TTa)Tr(T^aj  (4.28) 


1 


where  again  a  =  a/2  -  80  and  f,  =  2f,pZ. 


| 
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Expressions  for  Bj  and  Bj  valid  to  order  e2  can  he  obtained  by  using 
(4.27)  and  (4.28)  in  the  expansion  (4.1)  in  the  sane  fashion  that  the  expres¬ 
sion  (4.12)  for  B0  is  obtained  from  (4.1)  using  u0  =  b0.  We  do  not  write  out 
the  corresponding  expressions  for  B^  and  EL  here  but  only  note  that  the  functions 
b.  and  b.  deduced  from  (4.2)  are  obtained  under  the  restriction  |z|  >  |z  |  >0 

j  j  v  1  1  o 

and  within  the  argument  range  (4.15).  In  other  words,  B_.  and  B^  can  not  be  evaluated 
at  the  point  z  =  0  using  the  expansion  in  e2  technique.  The  problem  of  applying 
(4.1)  at  the  origin  for  the  solutions  B^  and  B^  arises  because  the  function 
U (a , ci  ;r )  in  the  expressions  for  b^  and  b^  is  singular  at  the  origin  (z  =  0)  for 

ct  >  1.  However,  the  behavior  of  the  solutions  B.  and  B  for  I  z  I  ->  0  for  all  values 

J  J 

of  a  can  be  found  by  employing  the  relations  (2 . 18) - (2 . 23) .  For  example,  multiplying 
(2.22)  by  ei2l,a-  and  using  (2.9),  (2.10)  and  (2.17)  yields 


A2(0,1)  =  -B^O.l)  +  el2Tfa  B3 


(4.29) 


Using  (2.23)  with  j  =  1  in  (4.29)  and  (2.19)  and  (2.22)  in  the  resulting 
expression  gives 

B3  =  [e12™  -  l]'1  (A2(0,1)  +  A:  +  A3  (0, 1)  -  B0}  .  (4.30) 

The  expression  (4.12)  can  be  used  to  evaluate  BQ ,  for  all  |z|  <  <«,  while  the  power 
series  expansion  (3.24)  can  be  used  to  evaluate  the  functions  A^  for  small  argu¬ 
ment  ( [ z |  <<  1).  Thus,  expressions  such  as  (4.30)  are  particularly  useful  in 
evaluating  and  EL  in  the  neighborhood  of  the  origin.  For  large  values  of  |z|, 
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on  the  other  hand,  and  are  more  conveniently  expressed  using  the  c2 
expansion  of  (4.1)  in  terms  of  and  b^. ,  respectively. 

Since  the  solutions  related  to  the  Kummer  function  U(a,a;f;)  contain  the 


combination  b^  -  b^  [c.f.  (4.21)]  it  is  useful  to  construct  the  combination 
B3  -  B3  in  a  form  that  can  be  evaluated  as  |z|  -*  0.  To  do  this  we  multiply 
(2.23)  by  e  ^2Tra-  and  use  (2.17)  to  obtain 


B3  =  B3 (0,-1)  -  Bo(0,-l)  .  (4.31) 

Multiplying  (4.30)  by  e  l2TTa-  and  using  (4.31)  in  the  resulting  relation  gives 

B3  =  [ei2Tr“  -  l]'1  (A2  +  Aj  (0,-1)  +  A 3  -  B0  (1,0) }  .  (4.32) 

Finally,  combining  (4.30)  and  (4.32)  yields  the  desired  combination 

B3  -  B3  =  [e127ra  -  l]'1  {(e12™-  -  1)[A2  +  A^O.-l)  +  A3  ]  -  (l  -  e12™+)B0) 

(4.33) 

Hxpressions  similar  to  (4.33)  involving  B2  and  B3  rather  than  B3  can  be  found 
by  using  (2. 18)- (2. 22) .  The  expression  (4.33)  can  be  evaluated  at  integer  values 
of  a  by  assuming  that  B3  -  B3  is  an  entire  function  of  a  and  then  using 
I, '  llospi  tal '  s  rule  to  obtain 

B3  -  B3  =  [  [A2  (0,1)  +  A3(0,1)  +  A3  +  B0 (1 , 0) ] / 2 

♦  (l/2ui)(el27ra-  -  l)[3aA2  +  DaAi(0,-l)  +  ^Ai] 

•»  (l/27ti)  (el2lTa+  -  iPaBol  I  •  (4.34) 


a=n 
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where  n  =  0,±1,±2,...  .  The  expressions  (4.33)  and  (4.34)  show  that  the  com¬ 
bination  B3  -  B3  is  finite  and  non-singular  at  the  origin  for  all  values  of  a. 

On  the  other  hand,  the  expression  (4.21)  with  j  =  3  indicates  that  b3  -  b3  is 
singular  at  the  origin  for  a  >  1.  Thus  the  exact  solutions  of  (1.3)  are 
rigorously  well  behaved  at  the  origin  although  the  restricted  expressions  obtained 
by  using  an  expansion  in  e2  indicate  otherwise,  which  is  a  consequence  of  the 
singular  perturbation  character  of  this  problem.  Expressions  for  B2  -  B2  and 
Bj  -  Bj  are  readily  obtained  from  (4.33)  using  (2. 18)- (2. 23) . 

Finally,  it  is  instructive  to  consider  a  special  case  in  which  the  parameters 
a  and  8  have  the  values 

a  =  v  ;  6  =  v/2 +  y  ,  (4.35) 

where  v  is  a  positive  integer  and  u  =  0,1,2,....  In  this  case 

a+  =  v  -  1  +  y  ; 

and 

a.  =  -(p  +  1)  . 

In  the  integral  representation  (2.7)  there  is  now  only  an  isolated  pole  of 
order  m  1  at  the  point  s  =  32*  The  points  s  =  ±g2  are  no  longer  branch  points 
and  the  s-plane  is  no  longer  cut  by  branch  lines,  so  that  B^(0,1)  =  B^  and  (2.23) 
then  gives 

B.  -  B.  =  Bn  .  (4.38) 

J  1  0 


(4.36) 

(4.37) 
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The  contour  B0  now  encloses  the  pole  at  s  =  g2  anc*  from  Cauchy's  Theorem  we 
immediately  have 


where  the  polynomials  are  given  by 


T,  -  e°»3/3  ^  <e't3/3> 

3  dtJ 


and  T0  =  1 .  The  expression  (4.43)  can  be  further  simplified  by  using  Rodrigues 

13 

Formula  for  the  Laguerre  Polynomials 


LV‘1+j  (x)  =  — _  eX  x1_V+j  — _  (xV_1+iJ  e'X) 

U-J  W  (U-J) !  6  X  dxp-j  (X  6  } 


Using  (4.45)  in  (4.43)  one  obtains  from  (4.39) 


B.  -  B.  -  B  *  2iri  (2g2)V_1  e"62*  e“a°3/3 
3  3  0  2 


f  p  (2o  )3 


for  p  =  0,1,2,...,  and  where  z  =  z  +  oq3/B 2.  Using  (4.44)  to  evaluate  terms 
in  (4.46)  results  in 


Bq  =  2iri  (2B2)V_1  e  &zZ  L^"1  (282z)  +  0(aQ3) 


Note  that  oq3  =  e2323»  so  that  the  correction  terms  to  Bq  as  given  by  (4.47) 


are  of  order  e2. 
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5.  Application  to  Physical  Problems 

In  applying  the  solutions  of  (1.3)  to  a  physical  problem  one  needs  to 

evaluate  the  solutions  A.,  B  ,  B.  and  B.  in  the  vicinity  of  the  real  axis 

J  0  J  J 

(Im  z  0+) .  However,  the  expressions  given  in  Secs.  3  and  4  for  some  of 
the  solutions  are  valid  only  in  certain  sectors  of  the  complex  z-plane.  Since 
some  of  these  sectors  contain  only  a  portion  of  the  real  axis,  the  expressions 
given  must  be  extended  by  means  of  analytic  continuation  to  include  the  entire 
region  of  physical  interest.  The  process  of  analytic  continuation  leads  to  a 
mixing  of  the  solutions  of  the  thermal  and  cold  plasma  classes.  This  mixing 
embodies  the  physical  phenomenon  of  mode  conversion. 

In  constructing  analytic  continuations  of  solutions  we  will  take  the 
imaginary  part  of  the  independent  variable  z  to  be  small  and  positive.  This 
assumption  leads  to  exponential  decay  of  waves  in  the  direction  of  propagation 
and  is  justified  on  physical  grounds  because  it  corresponds  to  adiabatic  switch-on 
of  the  exciter  at  frequency  w.  The  analytic  continuation  of  the  solutions  of 
(1.3)  can  be  accomplished  by  using  the  relations  (2. 18)- (2. 23)  which  are  derived 
by  applying  Cauchy's  Theorem  to  the  contours  defining  the  various  solutions. 

The  asymptotic  expressions  for  the  solutions  A^  as  given  by  (3.7)  are 
valid  in  sectors  of  aperture  4tt / 3  as  given  by  (3.4).  Since  die  solution  Aj  is 
defined  in  the  sector  5tt  /  3  <  8  <  3n  (z  =  |z|e  ),  which  includes  the  entire  region 
of  physical  interest,  it  needs  no  analytic  continuation.  Just  above  the  negative 


real  axis  0  =  3n  -  A,  where  A  is  a  positive  infinitesimal,  the  term  exp(-ic)  in 
the  asymptotic  expression  (3.7)  for  A;  varies  as  exp | C | .  Thus  A]  is  exponent i all\ 


growing  near  the  negative  real  axis.  Along  the  positive  real  axis  where  =  2n  +  A 
the  tern,  exp(-ic)  varies  as  exp  1  |  r,  |  so  that  A]  represents  an  outward  propagating 


wave  in  the  WKB  sense.  The  solution  A2  is  defined  in  the  region  it/3  <  0  <  5tt / 3 
which  includes  the  region  of  physical  interest  only  near  the  negative  real  axis 


where  0  =  x  -  A.  The  term  exp(-iO  in  (3.4)  for  j=2  then  varies  as  exp(-|^|)  so 
that  A 2  is  exponentially  decaying  for  Re(z)  <  0.  To  analytically  continue  A2 
we  use  the  relations  (2.18),  (2.20)  and  (2.22)  to  obtain  (see  Fig.  4) 


A2  =  -A3  -  B2  +  B2(1,0)  -  Ai(l,0), 


for  Re(z)>  0.  All  of  the  functions  on  the  right-hand  side  of  (5.1)  are  defined 
in  sectors  that  include  the  region  just  above  the  positive  real  axis  as  can 
be  verified  by  (3.4)  and  (4.15).  The  relation  (5.1)  indicates  that  the  solution 
A 2  undergoes  mode  conversion  when  the  Re(z)  changes  sign.  In  this  paper 
we  use  the  term  mode  conversion  to  denote  a  process  wherein  a  solution 
of  one  class  produces  a  solution  of  the  other  class.  For  example,  as  indicated 
by  (5.1),  a  solution  of  the  thermal  class,  A2,  near  the  negative  real  axis, 
leads  to  solutions  of  the  cold  mode  class  as  well  as  thermal  mode  class 
solutions  near  the  positive  real  axis.  Finally,  the  domain  in  which 
the  solution  A3  is  defined  includes  the  region  above  the  positive  real  axis  but 
not  the  region  above  the  negative  real  axis.  Near  the  positive  real  axis  0  =  A 
for  j=3  and  A3  varies  as  e  thus  representing,  in  the  WKB  sense,  an  inward 

propagating  wave.  Using  (2.18),  (2.21)  and  (2.20)  to  analytically  continue  A3 
we  find  that  (see  Fig.  4) 


a3  =  -Al (0,-1)  -  B3  +  B3 ( 1 ,0)  -  a2  , 
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where  all  of  the  functions  on  the  right-hand  side  of  t.5.2)  are  defined 
in  the  region  just  above  the  negative  real  axis. 

The  expressions  derived  for  the  solutions  and  are  also  valid  <mlv 
in  certain  sectors  of  the  z-plane  as  indicated  by  (4.15).  The  solutions  Bj  and 
B3  are  defined  for  the  entire  region  above  the  real  axis  and  need  no  analytic 
continuation.  The  solutions  B2  and  B2  are  defined  only  in  the  region  above  the 
positive  real  axis  and  must  be  continued  by  using  (2.21)  and  (2.22)  from  which 
we  obtain 


(5.3a) 

(5.3b) 

Both  functions  on  the  right-hand  side  of  (5.3a)  and  (5.3b)  are  defined  in  the 
region  just  above  the  negative  real  axis.  B2  and  B2  both  undergo  mode  conversion 
when  Re(z)  changes  sign  as  indicated  by  (5.3a)  and  (5.3b)  because 
they  generate  a  solution  of  the  thermal  mode  class  in  addition  to  a  solution 
of  the  cold  mode  class.  The  solutions  Bj  and  Bj  are  defined  for  the  region  just 
above  the  negative  real  axis  and  can  be  continued  to  include  the  region  just 
above  the  positive  real  axis  by  using  (2.18)  and  (2.19)  to  obtain 


B2  -  B3  +  Aj (0,-1) 

and 

B2  =  B3  +  Aj 


Bt  =  B2  +  A3  ,  (5.4a  1 

and 

Bx  =  B2  +  A3(0,1)  .  (5.4b) 


Again  the  solutions  Bj  and  undergo  mode  conversion  when  Re(z)  changes 
sign. 

In  addition  to  having  expressions  for  the  solutions  in  the  region  above 
the  real  axis  a  useful  set  of  four  linearly  independent  solutions  must  be 
selected  from  the  general  set  in  order  to  solve  a  well  posed  problem  of 
physical  interest.  In  particular,  we  have  in  mind  applying  the  solutions  of 
(1.3)  to  the  inhomogeneous  problem  in  which  a  driving  source  is  present  in 
the  plasma.  We  thus  shall  discuss  four  linearly  independent  solutions  that 
are  convenient  in  constructing  a  Green’s  function  for  (1.3).  The  linear 
independence  of  a  solution  set  can  be  formally  established  by  evaluating  the 
system  Wronskian.  Rather  than  perform  this  calculation  we  defer  it  to  a 
later  paper  and  only  present  here  a  heuristic  argument  for  linear  independence. 

A  Green's  function  must  satisfy  certain  boundary  conditions  and  we  impose 
these  beforehand  in  choosing  our  solution  set.  The  boundary  conditions  we 
impose  are  :  1)  wave- like  solutions  must  correspond  to  transport  of  energy  away 

from  the  source,  and  2)  solutions  must  be  bounded  as  |Re(z)|  ->  °°.  Since  the 
thermal  mode  and  cold  plasma  mode  classes  of  solutions  have  vastly  different 
scale  sizes  away  from  plasma  resonance,  our  strategy  for  identifying  a  linearly 
independent  set  is  to  select  one  pair  of  linearly  independent  solutions  from 
each  class.  We  are  then  assured  of  the  linear  independence  of  pairs  of  solutions 
from  different  classes  because  of  their  different  scale  sizes.  From  the  cold 
mode  class  of  solutions  we  choose  the  pair  B2,  B2  while  from  the  thermal  class 
of  solutions  we  choose  the  pair  Aj,  A2. 
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From  the  solution  pair  Bo,  Bo  we  next  construct  a  pair  of  linearly  inde¬ 
pendent  solutions,  one  of  which  is  bounded  as  z  ->-  °°,  while  the  other  is  bounded 
for  z  -*•  This  solution  pair  is 


BR  =  [(2b2)1-a/2TTi]  !'(-«_)  e  11Ta-{B2(a,31,62,e2,0 


B2(a,3j  ,B2,c2,.f:) } 


(5.5; 


B  =  [(2B2)1_a/2TTi]r(-a+)  e  1TO+{B? (a,-Bi , 3? ,-e2,-C)  -  B:,  (a  ,-P  j  ,  ;>  ,-  )■• 

L  ( 3 .  h  ) 

Note  that  in  (5.6)  we  have  used  the  transformation  (2.12)  so  that  we  are  assured 
Bj  is  a  solution  to  (1.3).  The  normalization  coefficients  in  (5.5)  and 
(5.6)  have  been  determined  using  (4.27)  and  (4.28)  and  are  chosen  so  that 
the  expressions  reduce,  when  e2  =  0,  to  the  simple  forms  given  below. 


B 


R 


e~U2 


U(a,a;F)  , 


(5.7) 


BI.  = 


+£/? 

e  ’  U (a  -  a,a ;-£) 


(5.8' 


where  a  =  </ 2  -  3  .  In  the  notation  of  Slater  these  two  solutions  arc 

o 

!  hen 


Br  =  e 


y;> 


(5.9) 


(5.10) 


so  that  using  (5.9)  and  (5.10)  together  with  (4.1)  we  find  that  the  IVronskian 

of  B  and  B  is 
K  L 


W(BR,BL)  =  e17r(0t/2+Bo)(r)  a  +  0 (e2)  , 


where  £  =  2B2Z-  In  (5.11)  we  have  again  assumed  that  Im(z)  is  positive  and  have 
ignored  contributions  to  the  Wronskian  arising  from  mode  conversion.  Mode  conver¬ 
sion  contributes  terms  to  the  Wronskian  that  are  a  product  of  cold  and  thermal  terms 
and  thus  are  rapidly  oscillating.  The  function  B  is  bounded  as  z  approaches 
positive  real  infinity  while  B  is  bounded  as  z  approaches  negative  real  infinity. 

L 

The  asymptotic  form  of  the  two  thermal  modes  Aj  and  A2  has  been  discussed 
before.  The  function  A2  is  exponentially  decaying  along  the  negative  real  axis 
while  Aj  is  exponentially  growing.  Thus  A2  and  Aj  are  clearly  linearly  independent. 
Near  the  positive  real  axis  \\  represents  an  outward  propagating  wave  while 
from  (5.1)  we  see  that  A2,  while  undergoing  mode  conversion,  contains  an  A3  term 
which  represents  an  inward  propagating  wave  so  that  Aj  and  A2  are  again  linearly 
independent . 

The  solution  set  B  ,  B  ,  A2 ,  Aj  is  convenient  for  constructing  a  Green's 

function.  If  we  denote  the  source  location  by  z'  then  the  solution  pair  B^,  A2 

can  be  used  exclusively  for  Re(z-z')  <  0  because  they  are  the  only  pair  bounded 

at  negative  real  infinity,  while  the  remaining  pair  BR,  A^  can  be  used  exclusively 

for  Re(z-z')  >  0  because  they  are  the  only  pair  bounded  at  positive  real  infinity. 

While  we  do  not  construct  a  Green's  function  here  we  can  use  the  solution  set 

B  ,  Aj ,  B  ,  A 2  to  illustrate  the  type  of  mode  conversion  that  may  occur  in  a 
R  L 


plasma  with  a  driving  source. 
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As  an  illustration  of  mode  conversion  in  a  physical  problem,  'onsider  the 

solution  B  for  Re(z)  >  0.  This  cold  mode  solution  as  given  by  (5.5)  depends 
R 

upon  B2  and  B2  which  can  be  analytically  continued  to  the  region  just  above  the 
negative  real  axis  by  using  (5.3a)  and  (5.3b)  from  which  we  obtain,  for  ile(z)  <  0, 


B  =  [  ( 2 32 ) 1  a/27ri]r(-a_)e  *  ~  (B3  (ct ,  Bj ,  B2  »£2  » £) 

R 


B3  (a ,  B  j ,  62  ,  t  ^  ,  t, ) } 


-  [  (2B2) 1  a/r(l  +  a_)]e"i27ra~  Aj  (a,3j  ,32,£2,0  . 


(5. 12) 


We  have  already  indicated  that  Aj  is  exponentially  decaying  near  the  negative 
real  axis  while  the  combination,  B3  -  B3,  is  defined  in  the  region  above  the 
negative  real  axis  and  according  to  (4.21)  is  proportional  to  the  hypergeometric 
function  U(a,cx;0  for  e2  =  0.  The  expression  (5.12)  illustrates  mode 
conversion  in  that  a  solution  of  the  cold  class  on  one  side  of  plasma  resonance 
consists  of  a  combination  of  cold  and  thermal  modes  on  the  other  side  of  resonance. 

We  note  that  in  this  particular  example  the  amount  of  mode  converted  thermal 
mode  is  proportional  to  1  / f ( 1  +  u_)  ,  and  thus  can  be  -zero  when 

1  +  u_  =  a/2  -  R  =  ,  (5.13) 


where  f  =  0,1,2,....  When  a  is  an  integer,  as  is  the  case  for  a  plasma, 
the  condition  (5.13)  is  identical  to  (4.35).  In  Section  4  we  found  that  under 
these  conditions  the  cold  mode  solutions  were  everywhere  proportional  to  B 
which  in  turn  could  be  expressed  as  a  sum  involving  Laguerre  Polynomial^  [see 
(3.4b)].  Thus  at  certain  special  values  of  the  parameter  R  the  phenomenon  of  mode 
conversion  does  not  occur  in  the  sense  that  the  amplitude  of  the  thermal  mode  is 
zero  and  furthermore,  at  these  values  of  R  the  cold  solutions  arc  proportional  to 
haguerre  Polynomials  for  c?  =  0. 
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The  quenching  of  the  mode  conversion  process  does  not  occur  for  all 
possible  solutions,  e.g. ,  the  solutions  B  for  Re(z)  <  0,  corresponding  to  a 
source  located  in  the  region  z  >  0.  Using  (5.6)  and  the  relations  (5.3a)  and 
(5.3b)  with  the  proper  parameter  values,  we  obtain  the  analytic  continuation, 
for  Re(z)  >  0, 

Bl  =  [(2f32)1_a/27ri]r(2-ta+)  e"1™*  (B3  (a ,-gj ,  g2  ,-e2  ,~0  -  B3  (a.-gj ,  g2  ,-e2  ,-?) } 

-  [(2t32)1_a/r(l+a+)]  e_l2lia+  Aj  (a ,-gj ,  62  ,-e2  ,-?)  .  (5.14) 

Since  a  and  6  are  positive,  l/r(l+a+)  =  l/r(a/2+B)  is  never  zero  and  mode  conver¬ 
sion  always  occurs  in  this  case.  Furthermore,  the  mode  conversion  process  is  not 
restricted  to  the  generation  of  thermal  modes  by  cold  modes.  The  thermal  modes 
also  undergo  the  mode  conversion  process  and  generate  cold  modes.  Although  the 
solution  Ai  does  not  undergo  mode  conversion  when  Re(z)  changes  sign,  the 
solution  A 2  does  as  is  illustrated  by  the  analytic  continuation  of  A2  given  in 


(5.1). 
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6.  Conclusions : 

By  applying  the  Laplace  integral  technique  we  have  obtained  integral 

representations  of  the  solutions  of  the  fourth  order  differential  equation  (1.3). 

The  solutions  are  distinguished  by  the  contours  along  which  the  integral  is 

evaluated.  The  solutions  A^  are  obtained  by  integrating  along  contours  having 

both  end-points  at  infinity  and  have  been  identified  as  belonging  to  a  class 

of  solutions  characterized  by  short  scale  lengths  and  thus  are  referred  to  as 

thermal  modes.  The  solutions  B  ,  B.  and  B.  are  obtained  bv  integrating 

°  .1  j  h 

along  contours  with  at  least  one  end-point  at  the  branch  point  s  =  -f ,  and 

belong  to  a  class  of  solutions  characterized  by  long  scale  lengths  and  are 

referred  to  as  cold  modes.  The  properties  of  the  exact  solutions  A  ,  B  ,  B 

J  o  j 

and  Bj  have  been  elucidated  by  expressing  them  in  terms  of  more  familiar 
functions  using  power  series  and  asymptotic  expansions. 

The  solutions  A^  have  been  evaluated  for  large  and  small  values  of  the  inde¬ 
pendent  variable  z  by  using  different  techniques.  A  power  series  expansion  in 

z  has  been  derived  in  (3.24)  in  order  to  determine  the  solutions  A  in  the  neighbor- 

J 

hood  of  the  origin.  The  solutions  A^  are  clearly  finite  and  well  behaved  in  the 
vicinity  of  the  origin.  For  large  values  of  z  the  functional  form  of  the  solu¬ 
tions  Aj  has  been  determined  in  (3.11)  using  the  saddle  point  method  of  integration 
to  obtain  the  leading  term  in  an  asymptotic  expansion.  For  real  values  of  the 
independent  variable  these  asymptotic  expressions  represent  either  exponentially 
decaying  or  growing  solutions  or  propagating  wave- like  solutions. 

Expressions  for  the  solutions  B0 ,  ,  and  have  been  obtained  by  expanding 

the  integral  representation  in  powers  of  e2.  This  technique  allows  an  expression 


to  be  obtained  for  the  exact  solutions  B0,  B^ ,  B^  in  terms  of  the  corresponding 
solutions  b0,  Ik,  Ik  of  the  second  order  differential  equation  obtained  from 
(1.3)  by  setting  c2  =  0.  The  expansion  (4.12)  obtained  in  this  manner  for  the 
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solution  B0  involves  the  Kummer  function  M(b,c;z)  and  is  valid  in  the  entire 
z-plane  (|z|  <  °°) .  In  particular  the  solution  BQ  can  be  evaluated  in  the 
neighborhood  of  the  origin  using  (4.12)  and  thus  is  useful  in  evaluating 
other  members  of  the  cold  mode  solution  class  near  the  origin. 

An  expansion  in  powers  of  e2  can  be  obtained  for  the  solutions  B^  and 
but  involves  the  restriction  |z|  >  |z0|  >  0.  The  expansion  involves  the 
corresponding  solutions  b^  and  b^  which  contain  the  Kummer  function 
U(b,c;z)  which  is  not  bounded  at  the  origin  for  Re(c)  >  1.  The  e2 
expansion  is  helpful  in  evaluating  the  solutions  and  for  large 
values  of  z.  The  solutions  B^  and  can  be  evaluated  in  the  neighborhood  of  the 
origin  by  using  Cauchy's  theorem  to  establish  relationships  among  the  various 
solutions  as  given  in  (2. 18)- (2. 23) .  Thus  the  solutions  B^  and  B^  can  be 
expressed  entirely  in  terms  of  various  combinations  of  the  solutions  and  B0 
as  shown  for  B3  in  (4.30).  Since  expressions  for  the  solutions  A^  and  B0 
valid  in  the  neighborhood  of  the  origin  have  been  obtained  the  solutions  B^  and 
Bj  can  also  be,  in  principle,  evaluated  there.  Unlike  the  associated  solutions  to 
the  second  order  differential  equation,  b^  and  ,  which  can  have  divergent  be¬ 
havior  at  the  origin  for  a  >  1 ,  the  solutions  B^  and  B^  are  always  finite  at 
the  origin.  The  physical  interpretation  of  this  result  is  that  the  inclusion  of 
thermal  effects  keeps  the  amplitude  of  the  electric  field  finite  at  plasma  reso¬ 
nance  through  the  production  of  short  scale  thermal  modes.  The  production  of  short 
scale  waves  near  plasma  resonance  involves  the  process  of  mode  conversion. 

For  certain  values  of  the  parameters,  a  and  3,  namely  a  a  positive  integer 
and  6  =  e2823/2  +  6 1 / 2 62  a  positive  half-integer,  the  topology  of  the  inte¬ 


gration  plane  is  greatly  simplified.  While  the  sector  structure  in  the  integration 
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plane  remains,  thereareno  branch  points  but  only  an  isolated  pole.  The 
quantities  ot+  and  a_  have  integer  values  under  these  conditions  so  that  the 
phase  factors  exp(i2irma+)  and  exp(i2ima_)  are  unity  and  as  indicated  by  (2.17) 
all  solutions  X(m,n)  are  equivalent  to  the  solution  X  evaluated  on  the  principal 
sheet.  The  combination  of  solutions  is  equivalent  to  the  solution 

Bq  as  indicated  in  (4.38).  Furthermore,  the  solution  B0  can  be  expressed  as 
a  series  of  generalized  Laguerre  Polynomials  as  given  in  (4.46)  which  is 
convenient  for  numerically  evaluating  quantities  for  physical  applications. 
Finally,  we  note  that  as  shown  in  (4.34)  the  combination  B3  -  B3  generally 
contains  an  admixture  of  the  thermal  mode  solutions  .  For  these  special 
parameter  values,  however,  the  combination  B3  -  B3  =  B0,  hence  it  contains  no 
thermal  mode  solutions.  Physically  this  indicates  that  mode  conversion  does  not 
occur  at  these  special  parameter  values.^ 

The  asymptotic  expressions  given  for  the  solutions  ,  B^  and  B^  are 
defined  in  certain  sectors  of  the  z-plane  as  given  by  (3.4)  and  (4.15).  Not  all 
of  these  sectors  contain  the  region  just  above  the  real  axis  which  is  the  region 
of  interest  in  physical  applications.  The  solutions  can  be  analytically 
continued,  however,  and  it  is  this  process  which  gives  rise  to  the  phenomenon  of 
mode  conversion.  A  solution  of  either  the  thermal  or  cold  mode  class  on  one 
side  of  plasma  resonance  gives  rise  to  a  combination  of  both  classes  on  the 
cither  side  of  resonance.  Physically  the  mode  conversion  process  serves  to 
limit  the  amplitude  of  the  solutions  at  plasma  resonance.  Specifically,  the 
solutions  of  the  purely  cold  plasma  (c2  =  0)  which  exhibit  singularities  at 
the  origin  now  correspond  to  solutions  which  are  finite  at  the  origin  but 
produce  thermal  modes  through  the  mode  conversion  process  which  carry  wave 
energy  away  form  the  resonance  region.  In  addition,  for  special  values  of  the 
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APPENDIX 


A. 1  General  expression  for  the  derivatives  of  u(z) 

The  nth  derivative  of  the  function  u(a,Bj ,62,£2»z)  with  respect  to  z, 
which  we  denote  as  u^n\  can  be  obtained  from  the  integral  representation  (2.7). 


u  n  (ay  ,32  »e2»z)  =  (-)11  /  sn(s-B2)°l~(s+B2)a't'  exp(-s3e2/3-sz)ds  , 

C 


(A. 1.1) 


where  C  is  any  of  the  contours  discussed  in  Section  2.  To  find  an  expression  for 
the  first  derivative  we  note  that 


s  =  (1/4B2)[(s  +  B2) 2  -  (s  -  B2)2], 


(A. 1.2) 


together  with 


u(a+2,6i+2B2 .B2*e2 »z)  = 


/  (s+B2)2(s+B2)a+(s-B2)a  exp(-s3e2/3-zs)ds  ; 
C 


(A. 1.3) 


u(a+2,Bi~2B2 .62 »e2 ,z)  =  /  (s-B2) 2 (s+B2)a+(s-B2)a“  exp(-s3e2/3-zs)ds  . 

C 


(A. 1.4) 


Using  (A. 1.1)  with  n  =  1,  the  expression  (A. 1.2)  together  with  (A. 1.3)  and  (A. 1.4) 


gives 


u^1\a,Bi,B2,e2,z)  =  (-1/4B2)  [u(a+2,  Bi+2B2  ,B2  ,e2  ,z)  -  u(a+2,h1-2;-.  ,>  '  ,z)  ]  . 


The  second  derivative  of  a  can  be  found  in  a  similar  fashion  by  using  the 


expression 


s2  =  (s  -  02)  (s  +  32 )  +  b22  ; 


(A. 1.6) 


we  then  obtain 


(2) 


(a,6l ,$2 ,£2 »z)  =  u(a+2,bi,B2,e2,z)  +  622  u(a,3j ,P2 ,z) 


(A. 1.7) 


Repeated  use  of  the  expression  (A. 1.7)  yields 


u(;n)(a,B1,62,e2,z)  =  £  f”|  B22  (n~j }  u(a+2j  ,  ,  B2 ,  z)  , 

j=0 


(A. 1.8) 


n  : 


where  n  =  0,1, 2, 3,...  and  " '  =  n!/(n-j)!j!,  are  the  binomial  coefficients. 

Application  of  (A. 1.5)  to  (A. 1.8)  then  yields 


n 


u(2n+1) (a,6i,62,e2,z)  =  (1/4)  [ 

j-o 


"  1  822(n_j)_1 


[u(a+2(j  +  l),Bi-262,82,e2,z)  -  u(a+2 (j+1) , B1+2B- ,32 ,r2  ,  z) ]  .  (A. 1.9) 


A . 7  A  series  expansion  of  u  in  powers  of  e2 

In  this  appendix  we  outline  the  procedure  for  obtaining  a  series  expansion 
:  ■  1 1  t . .  i  , ,l2,z)  in  powers  of  e2  and  explicitly  calculate  the  first  two  co¬ 
efficients  of  this  series.  We  first  assume  that  the  integrand  of  (2.7)  is  an 
m.ilvtic  function  of  rz  in  the  vicinity  of  e2  =  0  and  expand  a  portion  of  the 
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integrand  in  a  Taylor's  series  about  the  point  e2  =  0.  We  will  discuss  the 
validity  of  this  assumption  later.  The  integrand  of  (2.7)  depends  on  e2  in¬ 
directly  through  the  parameters  a+,a_  and  directly  as  exp[-e2(s3/3  +  |322s)]. 
Thus  part  of  the  integrand  in  (2.7)  can  be  written  as 


(s  +  g2)a+(s  -  e2)a'  exp[-e2 (s3/3  +  B22s)] 


(A. 2.1) 

*  (s  +  62)a+  (s  -  $2)a-  [l  -  e2{s3/3  +  B22s-  (823/2)£n['(s+82)/(s-82)]}j  + - 


where  we  have  explicitly  shown  only  the  first  two  terms  of  the  Taylor  series 
expansion.  The  parameters  a+  and  a  are  simply  the  values  of  a+  and  a_  at  e2  =  0, 
namely. 


a+  =  a+(e2  =  0)  =  a/2  -  1  +  Bi/2B2  J  (A. 2. 2) 

a  =  a  (e2  =  0)  =  a/2  -  1  -  Bi/2B2  ♦  (A. 2. 3) 


Using  (A. 2.1)  in  (2.7)  we  obtain 


u(a, Bi , 82 ,e2 , z)  =  /  ( s  -  B2)a-(s  +  B2)3+  e  ZS  ds 
C 

-  e2 (B22/  s(s  -  B2)a-(s  +  B2)a+  e‘ZSds  +(1/3)/  s3(s  -  B2)a*(s  +  B2)a+  e"zsds 
C  C 


-(S23/2)/»n[(s+(i?)/(s-p.2)](s+  B2)a+(s  -  B2)a"  e_ZSds)  +  .... 
C 


(A. 2. 4) 
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Noting  that 


/ ln[  (s  +  S2)/(s  -  S2)](s  +  62)  +(s  -  02 )  •  e  ds 
C 


=  2S2  ~  /  (s  +  02)a+(s  -  02)a“  e_ZS  ds  , 


(A. 2. 5) 


and  using  (A. 1.1)  together  with  (A. 1.9)  we  obtain 


u(a,6} ,$2,e2 ,z)  =  u(a,6i ,B2,0,z) 


-(e2/4B2)  ,[l/3]u(a  +  4,3i  +  2B2,S2,0,z)  +  C 1/3  +  322]u(a  +  2, Si  -  202,P2,O,z) 


-  [l/3]u(a  +  4,0!  -  2S2,62,0,z)  -  [l/3  +  322]u(ct  +  2, Si  +  232,82,0,z) 


4S25  JL  u  (a 


,61,62,0,z)J  +  ... 


(A. 2.6) 


Writing 


u  _  u (a , Si , S2 , 0 , z)  =  /  (s  -  Bo) 3  (s  +  S2)  +  e  S"  ds  , 

c 


(\.2. -) 


we  can  then  express  (A. 2.6)  as 


u(«, 0! , S2 ,z)  =  u(a,hj ,S2,0,z)  +  t  '  (  1  u) 


=  uq  -  e'-hlu  )/4B2  +  ...» 


(  A  .  2 . 8  ) 


-52- 


where  h(u  )  is  a  functional  of  uq  and  denotes  the  term  in  braces  in  (A. 2. 6). 

That  is,  once  the  functions  u^  (which  are  independent  of  the  magnitude  of  t) 
are  known,  the  behavior  of  u(a,  Bi , 62  »e""  > z)  can,  at  least  in  principle,  be 
determined  for  any  value  of  e2  provided  the  Taylor  series  expansion  is  valid. 

For  the  Taylor  series  to  be  valid  the  functions  uq  must  be  analytic  functions 
of  e2  in  the  neighborhood  of  e2  =  0.  Furthermore,  the  coefficients  in  the  Taylor 
series  must  be  bounded  for  all  values  of  z.  The  functions  uq  are  analytic 
functions  of  e2  only  for  the  contours  corresponding  to  the  solutions  Bq ,  and  Bj 
because  these  solutions  become  solutions  of  the  second  order  differential  equation 
when  c2  =  0.  On  the  other  hand,  u  integrated  along  the  contours  for  Bq ,  B  and 
B  can  be  singular  for  certain  values  of  z  and  the  coefficients  of  the  Taylor 

j 

series  are  not  bounded  as  z  approaches  these  points.  In  the  problem  treated  here, 
for  example,  singularities  in  b  and  b  can  occur  at  the  point  z  =  0  for  a  >  1. 

In  this  case  the  expansion  in  powers  of  ez  is  valid  only  for  z  bounded  away  from 
the  origin,  i.e.,  |z|  >  jzj  >  0,  for  a  given  zq,  and  for  arg  z  restricted  to  the 

appropriate  sectors,  as  given  in  the  text  by  (4.15). 


Figure  1.  Contours  corresponding  to  solutions  of  (1.3)  must  proceed  to  infinity 
in  the  shaded  open  sectors  labeled  1,  2  or  3.  The  integration  plane  is  cut  along 
the  real  axis  from  6 2  t0  positive  infinity  and  from  -fs2  to  negative  infinity.  The 
plane  of  the  paper  represents  the  principal  sheet  with  m  =  0,  n  =  0  in  (2.16). 

Figure  2.  The  contours  for  the  solutions  Aj ,  A2  and  A3  begin  and  end  at  infinity. 

The  contour  for  the  solution  A2  passes  onto  the  adjacent  Riemann  sheet  with  m  =  1, 
n  =  0  in  (2.16)  when  crossing  the  branch  cut. 

Figure  3(a).  The  contours  for  the  solutions  B0  and  have  at  least  one  end 
point  at  s  =  -£S2.  Contours  that  cross  the  branch  cut  pass  onto  the  adjacent 
Riemann  sheet  on  which  m  =  0  and  n  =  1  in  (2.16)  and  are  shown  dashed.  The  contour 

for  the  solution  B3(0,1)  starts  on  the  sheet  (0,1)  and  crosses  onto  the  principal  sheet. 

Figure  3(b).  The  same  as  Fig.  3(a)  but  for  the  contours  corresponding  to  the 

solutions  B . . 

J 

Figure  4.  Contours  for  all  of  the  solutions  shown  in  Figs.  2  and  3  are  combined 
to  facilate  the  derivation  of  the  relations  (2. 18)— (2 . 23) . 

Figure  5.  The  position  of  the  saddle  point  used  to  obtain  an  asymptotic  represen¬ 
tation  for  the  solutions  A^  is  shown  as  a  function  of  0  =  arg(n).  The  saddle  point 
starts  just  above  the  right  hand  branch  cut  for  0  =  -it  and  traces  out  a  circular 
path,  crossing  onto  the  adjacent  m  =  1,  n  =  0  Riemann  sheet  to  end  just  below  the 

branch  cut  at  0  =  3tt.  The  contour  for  the  solution  Aj  lies  in  the  (1,0)  Riemann 


sheet . 
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